RESONANT MODE INTERACTIONS 
IN NONLINEAR LATTICES 



V.V. Konotop 
Departamento de Fisica and CFMC, 
Universidade de Lisboa, Complexo Interdisciplinar, 
Av. Prof. Gama Pinto 2, Lisbon 1649-003, 
Portugal 

February 1, 2008 

Abstract 

Resonant mode interactions in weakly nonlinear multi-dimensional lat- 
tices and related effects are described. We concentrate on formal descrip- 
tion of the phenomenon and consider as examples mode interactions and 
evolution equations for quadratic solitons in two-dimensional lattices, and 
coherent structures in a double-chain as examples of the application of the 
theory. 

1 Introduction 

Elastic interactions among quasi-particles, like phonons, electrons, magnons, 
etc., in solids constitute fundamental phenomena allowing one to observe mi- 
croscopic effects on the macroscopic level. Today they are very well studied 
and are described in numerous text-books (see e.g. Modelung, 1981). Although 
interactions mean nonlinearity and thus possibility of a complex dynamics such 
phenomena have been discussed, so far, on the language of quasi-particles. This 
means that usually only two-body interactions are taken into account. Mean- 
time, due to fundamental discoveries of sixties and early seventies, namely due 
to construction of the modern soliton theory, it became clear that the dynamic 
of nonlinear systems is much richer than one provided by the theory of two-body 
interactions. Namely, it became clear that together with linear quasi-particles 
one should consider solitons (or solitary waves) which are essentially nonlinear 
objects and cannot be constructed from the linear theory by any perturbation 
technique. In particular, speaking about nonlinear lattices one should consider 
not only phonons and interactions among them, but also solitary waves, called 
envelope solitons (Tsurui, 1972). Meantime solitons are approximate solutions 
of the nonlinear lattices and they are necessarily accompanied by generation 
of higher harmonics (Jimenez, 1999). The fact that corresponding dynamics 



1 



displays indeed a complicate behavior has been reported in (Economou, at al 
1994; Kopidakis, at al. 1994). The effect of high harmonics becomes especially 
important when the resonant conditions 

W3 = wi+a;2, q3 = qi + q2 + Q (i) 

or 

cj2 = 2cji, q2 = 2qi + Q (2) 

are satisfied. On the language of solid state physics one can interpret the first 
and second equalities in ([^) or in as conservation laws of energy and momen- 
tum for phonons, respectively. Using the terminology of the nonlinear optics 
these formulas can be called matching conditions or conditions of resonant three 
mode generation and resonant generation of the second harmonic. The present 
work describes various phenomena in nonlinear multi-dimensional lattices when 
one of the above resonance conditions is satisfied. 

The theory we develop is applied to lattices with a complex cell, to vec- 
tor lattices (with a simple or complex cell), or to couple lattices of different 
physical nature (like for example to a system describing electron-phonon inter- 
actions (Konotop, 1997)). Also, subject to certain modification the theory can 
be adjusted to describe coupled nonlinear chains. In the context of the present 
work the last system is interesting since it does not belong either class of one- 
dimensional nor two-dimensional lattices. Also a double chain is an example 
of a system when one can provide resonant conditions between localized and 
extended degrees of freedom. 



2 Properties of a linear lattice. 

Let us recover some general properties of a linear A/'D lattice [hereafter M {= 
1, 2, 3) indicates the dimension of the space]. It will be assumed that the lattice 
is finite with Nj (j = 1, ...,J\f) cells in jth direction, so that the total number of 
cells is A^ = A^i • • • N^/. It also will be assumed that each cell has a generalized 
degrees of freedom: they can be related either to different atoms, or to different 
components of the displacement vector. After all we will be interested in the 
limit Nj — !■ cx). We start with a generic case, when the evolution equation is 
written as follows 

*(n) = -5]j(n,l)*(l). (3) 
1 

Here 5* (n) is a.J\fxJ\f matrix with A/" = aJ\f having elements ipij (n) , n is an J\f- 
dimensional vector indicating the position of a cell: it can be represented in the 
standard form n = ''^j^jj where arc lattice vectors with rij = 0, 1 . . . Nj — 1, 
J(n, 1) is a A/" X A/" matrix of force constants which is considered to be hermitian, 
Jij(n, 1) = Jji(ji, 1) (the overbar stands for the complex conjugation) and having 
only positive eigenvalues, and a dot standing for the derivative with respect 
to time. It will be assumed that the second index in ipij refers to a branch 
of the spectrum, while the first one is associated with a generalized degree 
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of freedom. Lattice (||) is considered subject to periodic boundary conditions 
*(n + NjBj) = ^'(n) for all j. 

Evolution equation (||) can be associated with the spectral problem 



$(n;q)A(q) -^J(n,l)ci>(l;q) 



(4) 



where $(n;q) = ($i(n;q), . . . ,$^(n;q)): $^(n;q) are column eigenfunctions, 
A(q) = diag[Ai(q), . . . , Aj^(q)] , is the eigenvalue matrix: a;^(q) = A^(q), (r/ = 
1 . . .Af) and q is a vector in the reciprocal space. It will be the assumption that 
in the spectrum of the linear eigenvalue problem there are no degenerate points 
i.e. A^(q) ^ A^'(q) for any qii-q^rf. 

We will be interested in translationally invariant lattices, this is J(n, 1) = 
J{\ n — 1 I). This property allows one to represent the eigenfunctions in the form 

$(n, q) = A(q)e*i" (5) 

where A(q) is a matrix depending on the wave vector only. It solves the equation 

A(q)A(q) - J(q)A(q) (6) 

with J(q) = X^n '^(l ^ |)e~*''". It will be assumed that the eigenfunctions 
(}ii^(n, q) make up an orthonormal complete set. Thus 

^<ft(„.q')$(„.q)^A(q,q'), ^ $(n; q)$t (1; q) ^ A(n, 1) (7) 



where A(q, q') = /5q,q', (5q,q' is the Kronecker delta, / is a unit 4x4 matrix and 
the dagger stands for the Hermitian conjugate. It follows from (0) that A{q) is 
a unitary matrix A''{q)A{q) = I which means that <I>(n,q) is a unitary matrix, 
as well. Now the dispersion relation of the linear lattice can be expressed as 



(8) 
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Next we express the group velocity through the eigenfunctions of the linear 
spectral problem. To this end we consider a small variation q — > q + dq and 
the resulting changes of the eigenfunctions 



<i>(l;q) -^$(l;q) 



$(l;q) 



This variation can be interpreted as one caused by the perturbation 



(5J(n,l) = 



z(l-n).dq--((I-n).rfqf 



J(n,l) 



of the Kernel J(n, 1). Thus on the one hand the Taylor expansion yields 
A(q + dq) = A(q) + (dqV,)A(q) + ^{dqV,)'A{q} + Oi\dq\') 



(9) 



(10) 
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where = 9/9q. On the other hand we obtain from the standard perturbation 
technique 

A(q + dq) = A(q) + A(l)(q) + A(2)(q) + 0(|dq|3) 
where the diagonal matrices A^^' (q) have elements as follows 

A^'.^q) ^^J2idq■in~ l))if,,(n, 1; q) (11) 

l,n 

= E(^q • - 1; q) 

l,n 

-EE - - 1')) E (12) 

11;/ ^ ■^n a 

l,n 1 ,n CK7^?7 

where we have introduced the matrix 

i^^„(n, 1; q) = <i>t (n; q) J(n, l)<i>„(l; q) (13) 
Comparing (|ll|) and (12) with (|^) we compute the group velocity 



Vr,(q) = VgCJ,,(q) = ^, , V (n2 - ni)i^,p,(ni, na; q). (14) 



ni,n2 

and the tensor of the group velocity dispersion 
1 duj„ dujn 1 



J? = - 

ujr, dqi dq-j uj, 



+ ~E("» - - ^j)^w(n,i;q) 



l,n 

' 1'-, 



- - E E("^ - '^)K - E (15) 

l.n l',n' a#r, ^" 

3 Nonlinear model 

The subject of our main interest is the resonant mode interactions in a nonlinear 
lattice possessing quadratic nonlinearity. As in the previous section we start 
with the general model whose dynamics is governed by the equation 

*(n) = - ^ J(n, ni)«'(ni) - ^ J3(n, m, n2) : ^'(ni)*(n2) (16) 

ni ni,n2 

Here 4'(n) =col('0i(n), . . . , ,7/1^(11)), J(n, ni) is a A/" x A/" matrix possessing 
the properties described in the previous section and J3(n, rii, n2) is understood 
to be a tensor with components J3 jjua ' such that the last term in the formal 
expression in (^^ means that the jth component of the nonlinear term is given 

by 

(j3(n,ni,n2) : ^'(ni)^'(n2)) = ^ J3JJ1J2 (n, m, n2)V'ii (n^V-js ("2) 
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Below we concentrate on the case of nearest neighbor interactions. 

In order to describe mode interactions we employ the multiscale analysis. 
To this end we introduce a formal small parameter e, e ^ 1, which generates a 
set of temporal variables ti, = e'^t, (z/ = 0, 1, . . .) and a set of spatial variables 
no and = e'^n, (i/ = 1, 2, . . .) which are regarded as independent, will be 
considered as a continuous variable. Respectively, we look for the solution of 
( p^ in the form of the expansion 

vl/(n) = e*(°)(no, {rj; {U}) + e'^d) (no, {rj; {U}) + ■■■ (17) 

where {t^} and {t^} stands for the whole set of variables. Being interested in 
the process of the three-mode interactions we represent 
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V'(°Hno,K};{t,}) = ^^„(ri;ii)</.„„(no;qa)e-'"-('i°)*° +C.C. (18) 



Here, 0,,^(no;qQ) are the eigenfunctions of the eigenvalues uj^^ linked by (|l|), 
Aa{vi] ti) can be interpreted as the amplitude of the envelope of the a-th mode 
and "c.c." stands for complex conjugation. In a generic situation the ampli- 
tude A depends on all "slow variables" (ri, r2, . . . ; <i, t2, • • ■)• However being 
interested in the process of three wave interaction which corresponds to the 
scales defined by ri and ti we indicate only these variables in the argument 
of the amplitude. The scaling imposed implies neglecting the group velocity 
dispersion. 

Substituting the new variables in evolution equation (^) and gathering 
terms of the same order of the parameter e one recovers the dispersion rela- 
tion (H) in the first order. In the second order with respect to e subject to the 
matching condition (^) one obtains the following system of equations for slowly 
varying amplitudes 

,V.)A=*^™^^i^A^3 (19) 



dti 



'Vl 



^ + (V,. . Vr)A, = ^^™iHilSE)^,^3 (20) 

Oil UJr)2 

Oil ' ' 



where = -j^- The nonlinear coefficient Xr;ij;2i)3 (lii ^2) is given by 
Xriiri2ri3 (qi,q2)= •^3jii2j3(ni,n2,n3) 

"1,12, "3 jl J2 J3 

X'^ji,m(ni'qi)^i2,»)2(n2,q2)0j3,'73(n3,qi +q2) (22) 

where 0j_^(n,q) are entries of the column eigenvector $,,(n,q), j = 1,2, ...jN, 
(hereafter vectors of the reciprocal lattice Q are dropped in the arguments of 
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the functions). In the particular case of the second-harmonic generation, when 
conditions (|^) are satisfied one has 

^ + (v,,V,.Mi = ^^A^2, |^ + K,V,M2-^^? (23) 

Uti ^rii Oil ^ri2 

where for Xnirn = Xnirni'^ one must use the formula 

Xmm^ E E "^3JiJiJ2("i'"2,n3)(<^j,,^,(ni,q))^</>j,,^,(n2,2q) (24) 

ni,n2,n3 ji ,j2 



4 Second harmonic generation in a two-dimensional 
lattice. 

As the simplest example of the above approach we consider a 2D diatomic lattice 
having atoms with masses Mi and A/2 with linear interactions between nearest- 
neighbors and subject to a nonlinear on-site potential. In order to introduce 
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Figure 1: The example of a diatomic lattice. Solid lines indicate inter-atomic 
interaction taken into account. 

atom numbering we can choose a cell as is shown in Fig. |l]. Then the lattice 
Hamiltonian reads 

+ ^[(ui(n) - U2(n))2 + (ui(n) - U2(n + ei))^ 

+ (ui(n) - U2(n + 62))^ + (ui(n) - U2(n + ei + 62))^] 

+ ^[ul,{n)+ul,{n)] + ^[uli{n)+uUn)]'j (25) 

Here K2 is a linear force constant, K2j and K^j are respectively the linear and 
nonlinear constants of a force applied to the j-th atom in a cell and originated 
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by the on-site potential, Uij is the z-th component of the displacement of the 
j-th atom. We assume that the lattice is square with a unitary lattice constant: 
BiCj = 5ij. Then in order to represent the evolution equation in the form ( p^ 
we define ip{n) = col {y/MiUn, V-^2Ui2, ^^-^2^22) and calculate the 

components of J2 and J3. Then the dispersion relation for the linear lattice can 
be written in the form 



1 + v/i-7^^(q) 



(26) 



where a — [^^] , here [•] stands for the integer part, u = 2 ^'2 ' ^ 
F(q) = sm2(|) + sin2(|) + sin2(|) sin2(|). 



(27) 



4m 



(1 - m)2(fc2 + 2ka) + (1 + m)2 ' 



7^ (l + m)(l + /c^) 



ujI^ = 2^^2(1 + + = ™ = ^ and = and = 1,2 

correspond to the low frequency branches which are reduced to the accoustic 
modes at K^j = 0, while ry = 3,4 denotes high frequency optical branches. 

Peculiarities of the second harmonic generation in the 2D case follow from 
the dispersion relation. Firstly, in contrast to the ID case the resonant condi- 
tions (El) are satisfied on a continuous line which can be seen from the fact that 



OwithX = sin^(f ), 



the equation for the line can be written as J^ij CijX^Y^ 
Y = sin2(^) and constants Cy which are trivially found from (|2^). Secondly, 
energy transfer can occur between modes with the same direction of velocities 
or with different directions of velocities. In particular, the velocities can be 
orthogonal. The former case in many aspects is analogous to the second har- 
monic generation in the ID chain (Konotop 1996, 1997). Finally, it is possible 
to have simultaneous processes when one low-frequency mode can resonantly 
interact with two other modes. In particular the processes 2wi(qi) = W2(2qi) 
and 2a;i(qi) = W3(2qi) can occur simultaneously for qi = (tt, 0) if the lattice 

fe?+2fci + l , , _ 12fe?+21fci+10 



parameters are chosen as following m = 12- 



and k2 = 



'3fei+6fci+4 '"-^ 3(l+fci) 

Although in all other points the equations ( P3| ) are similar to the equations 
describing the second-harmonic generation in nonlinear optics, at such points 
equations (|2^ ) are not valid any more. 

The effective nonlinearity controlling SHG is now given by 




X3lA;,i(q)A„4(2q) + X32A;,3(q)A„3(2q)] 
i^3l<,2(q)A,2,2(2q) + X32<,4(q)A2.4(2q)] 



(28) 



Here yl^j (q) is the amplitude which can be trivially found from 
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If v^^ is parallel to v,,^ then this reduces the system ( p3| ) describing the 
evolution of the envelope to an effectively one dimensional problem 

= ^+-.^=^^-4? (29) 

Oti ox Uli Oti ox UJ2 

where x is the slow spatial coordinate in the direction of the group velocity 
vector and x — Xj)i))2(i)- This system is well studied in the context of nonlinear 
optics (see e.g. Suhorukov. 1988). It also reduces the 2D problem to quasi-one- 
dimensional one (Konotop, 1996). 

4.1 Case of orthogonal velocities 

If the velocities of the first and second harmonics are orthogonal then system 
( p3| ) is reduced to 

A1A2, -7T—+V2^-=l A{ (30) 

oti ox LUl oti Oy L02 

where x and y are spatial coordinates in the direction of the group velocities of 
the first and second harmonic, respectively. 

In a special case the system ( |30| ) can be linearized, using the ideas of (Suho- 
rukov. 1988). Indeed, let us introduce normalized amplitudes 

ai = -^e~'"Mi, 02 = ^e"*"M2, 2ai - aa = 7r/2 
V2wi <~oi 

and look for a solution of ( ^0|) such that the functions oi and 02 are real. Let 
us also define operators 

L ^ 9 L ^ ^ 

^ dt ^ dx^ ^ dt ^ dy 

Then system (^0[) can be rewritten in the form 

I/ifli = 0102, -^202 = -a^. (31) 

Next we apply Li to the second of the above equations. As a result we arrive 
at £2(^1^2 — 0-2) = 0. This means that Lifla — — C{x,y — V2t) where C(-, •) 
is a function which is determined from the boundary (initial) conditions. Then 
the linear equation for 1/ai takes the form 

Ll— + C{x,y-V2t)—^0. (32) 
ai ai 

4.1.1 Solitary wave solution 

The first simplification of (^) can be achieved by assuming that C(-, •) does not 
depend on x: C(-, •) = C(^2)- Here we introduce running variables: £_i — x/vi — t 
and ^2 = y/v2 — t. Then one can look for the solutions depending on the running 
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variables only: aj = aj i , C2) • Here we represent one of them which is generated 
by the Lorentz pulse C{^2) =const/(l + C2/''^2)^- It reads 



X\/^ + "^l [^i'^l^^i sinh(arctanf2) — coshfi cosh(arctanf2)] 
^ iiui 1^2 cosh(^i — arctan^2) ~ sinh(^i — arctani^2) 

^ X (1 + Cl) cosh(^i - arctan^2) 

The first harmonic represents a solitary wave pulse which is exponentially 
localized in the direction of the group velocity of the first mode and displays 
power-law-localization in the direction of propagation of the second harmonic. 
The second harmonic is dclocalizcd in the direction of the group velocity of the 
first mode. In that direction it has a kink-like shape. At t > the two modes 
being coupled with each other move without distortion in the direction defined 
by x/vi + y/v2 = 0. 

4.1.2 Spatial structure 

In order to illustrate the possibility of existence of stationary spatial structures 
we consider solutions of [or (|3l|)] in a space defined by a; > 0, y > when 
— subject to proper boundary conditions and subject to the assumption 
that V1V2 > 0. It reads {xq and yo are real constants) It reads 



^2ix,y) 



X {e^+^° sinh(?/ -I- yo) + e~^~^° cosh(j/ -I- yo)) 
viUJi (e~^~^o cosh(y + yo) ~ e^+^" sinh(j/ + j/q)) 
X (e^+^f sinh(j/ + yo) + e^^^^" cosh(y + yo)) 



The respective solution represents energy transfer form the first harmonic, lo- 
calized in the direction of the group velocity of the first mode, to the second 
harmonic, in the same direction having a kink-like shape, which occurs in the 
direction of the group velocity of the second mode. 



5 Quadratic solitons in a lattice 

Consideration provided in the previous section does not take into account the 
dispersion which is intrinsic property of the nonlinear lattices. This was justified 
by the choice of the scaling of the problem. Let us now turn to the case when 
the group velocity dispersion appears to be an important factor (Konotop and 
Malomed, 2000). Namely, assuming that resonant condition (^) is satisfied let 
us look for the solution of nonlinear model (|l^) in the form of the expansion 
[c.f. 0] 

*(n) = £2.1,(0) (no, {r,}; {U}) + e^vj/d) (no, {r J; {U}) + ■■■ (33) 
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Then the substitution 

2 

^(")(no,K};{t,}) = ^^„(ri;ii)</.„„(no;qa)e-'"-('i°)*° +C.C. (34) 

results in the evolution equations for the amplitudes (they are obtained in the 
third order of the asymptotic expansion) 

Some particular solutions of these equations (however in applications to the 
nonlinear optics) have been considered in (Malomed et al. 1997). 



6 Dynamics of a double chains 

In the present section we deal with phenomena originated by resonant interac- 
tions of longitudinal and transverse degrees of freedom in a system consisting 
of two coupled vector atomic chain depicted in Fig. |^. 




M-chain 



D-chain 



Figure 2: Schematic representation of a double chain. The axes x and y will be 
chosen along a and b, respectively. 

The potential energy of the interaction between the neighbor sites, say n 
and n + 1, is considered to he Ul{\sl + u„_|_i — u„|) in the u-chain and C/L(|a + 
v„+i — v„|) in the w-chain. The nearest neighbors in the both chains interact 
with the energy C/rdb + u„ — v„|). 

Thus the hamiltonian of the system can be written down in the form 

^-^ 2m 2m 

n \, 

+ C/L(|a + v„+i-v„|) + C/T(|b + u„-v„|)} (37) 
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The vectors u„ — {ux,n,Uy.ii) and v„ = ivx,n,Vy,n) are displacements of the 
atoms from their equihbrium positions in the u- and w-lattices, respectively. 
Without restriction of generality, in what follows the atomic mass is considered 
to be one: m = 1. 

If uju and iOy are a characteristic frequencies of excitations propagating along 
u- and v-chains. and LVt is a frequency of the transverse oscillations, it is natural 
to expect enhancement of the interaction when ntUt + + nycu^ = 0, nt, 

Uu, and Uy, being integers. The "lowest" resonant phenomena will be related 
to resonant three wave interactions either ujt = 2lui or toi = 2LUt. 

We will concentrate on excitation of a relatively small amplitude. Then the 
potential energy in (|37|) can be expanded in the Taylor series 



C/i(|a + u„+i-u„|) = l-K2{Ly{n)~af + lK3{Ly{n)-af + ... 

2 b 

C/T(|b + u„-v„|) = ig2(L(„)-6)2 + ig3(L(„)_5)3 + ... 

where 

1 /2 

Luin) [{a + Ux,n - Ux,n-lf + [Uy^n - U.y^n-lf] 

is the distance between two neighbor sites in the u-chain (L„ is defined by 
analogy) 

1 /2 

L{n) = [{Ux^n - Vx^nf + + "y,™ - 

is the distance between neighbor sites u- and v-lattices and the force coefficients 
are given by Kn = ^-g^r^, and Qn = ^j^J-''^ (if 2 and Q2 being positive). 

Since we are dealing with the phenomena originated by resonant three wave 
interactions the terms due to the nonlinearity of the central-force potentials Ul 
and Ut are not important in the leading order, since the lowest order nonlinear- 
ity results from the expansion of Lu{n), Ly{n) and L{n) with respect to small 
displacements of atoms from the equilibrium positions. Then the equations of 
motion read 



Q2, 
b 

-^[{Uy,n+l - Uy^n f - {Uy^n-1 - Uy,n)'^] (38) 



+ ^[{Vy,n+l - Vy^n)'^ - {Vy^n-l - Vy^n)"^] (39) 

K2 

[(^a;,n+l '^x,n')i^'Uy,n-\-l ~^ '^y.n) ~ ~ ^a;,n — 1 ) (^y,n ^y,n— l)] 

(40) 



11 



(41) 



6.1 Resonance cut — 2uji. 

As an example consider the situation when motion of u- and v- chains is syn- 
chronized, i.e. and respectively Wt = 2uji. As far as lattices are 
considered identical the above assumption means that the wave vectors of the 
carrier wave (cw) of excitations propagating along both chains coincide. We 
will concentrate on the case when the wave vector of a single lattice borders the 
boundary of the Brillouin zone. Then solutions of the hamiltonian equations 
for the lattice displacements can be represented in the form 







4^)(ti;xi)(-ire- 
















-e^ Ui)(ii;ari)e-2''^*' 



c.c. 



U,2 



+ C.C. 



Here we introduce a hierarchy of slow times tj = eH and spatial scales rij = 
(in {j = 0,1,...). The scales nj {j > 1) are slow ones and that is why we 
introduce Xj = arij (j > 1) which are regarded as continuous variables. The 
symbol O(e^) stands for all the terms of order of which have frequencies either 
different from oj for longitudinal components or different from 2uj for transverse 
component. Also in the expansion wc use the convention that in the arguments 
of the functions only the most "rapid" scales are indicated [so, for example 
Ui(ti]Xi) = ui{ti, 1,2, xi, X2, ...)]. "c.c." stands for the complex conjugation, 
i.e. Quantity 2w{^;t2, X2) can be interpreted as an is an enlargement of the 
distance between neighbor atoms in different chains which oscillate with opposite 
phases. 

In order to get solitonic effects this expansion should be provided up to the 
terms of order of e^. After that the expansion for u, v must be substituted in 
the Hamiltonian equations and terms of the same order of the small parameter 
must be collected. Dropping this straightforward standard calculus we pass to 
the description of the equations governing the dynamics of the double chain. 

In the lowest order with respect to the small parameter (i.e. in the order) 
the chains appear to be independent. The respective terms of the expansion 
of the longitudinal and transverse components give the two values for the fre- 
quency u) = 2\fK^ and w = \/2Q2, which must coincide. In other words the 



phenomenon can be observed subject to the condition 



da 



Sub- 



ject to the above conditions, the equations of the third order with respect to 
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epsilon become zero identically. Thus, computing the terms of order of one 
arrives at the system of equations 

dU d^U - - dV d^V 

-f^ + ^(C/-^f-0 (43) 

where U{t,xi) = ^ui^\t2,xi), V{t,xi) = ^vi^\t2,xi), W{t,xi) = ^w(t2,Xi) 
and T = y/K^t2/2. 

Let us consider a particular solution of the above system, which corresponds 
to out-of-phase oscillations of the atoms in the u- and v-chains: U = —V. Then 
the system (|2|)-(|^ takes the form 

dV d^U . .dW 



-i^ + ^+^UW = 0, -t— + 2U^ = (44) 

OT Oxf OT 

We notice that the hamiltonian governing this system is as follows 

^ = J (iC^o^P + 2I7V + 2WC/2) dx (45) 
System ( ji^ ) admits solutions in a form of two coupled dark solitons 

C/ = /32e2*/3'ttanh(/3x), W ^ -^e^'^'Ha.nh'ipx) (46) 
and two coupled bright solitons 

U = - - W^-- - (47) 

V2cosh(/3a;)' 2 cosh2(/3x) ^ ^ 

These solutions describe localized in space narrowing and broadening of the 
double chain in transverse direction. The respective trajectories of the atoms 
are parabolas which arc modulate by the envelopes (Eq) and (E^) respectively. 



7 Conclusion 

In the present paper we have considered several effects in multidimensional 
lattices which are originated by resonant wave interactions. All the examples 
were given for nonlinear Klein- Gordon like lattices. The theory however allows 
for straightforward modification for the sake of description of lattices of the 
nonhnear Schrodinger type (see e.g. Konotop 1997; Konotop, et al. 1999). 
Also, resonant muti-wave interactions in lattice can be considered in the similar 
way. 

Other relevant effects we did not discussed here are associated with the phase 
mismatch (i.e. when, for example, (|l|) is substituted by c^a = uji + uj2 + Acj, 
Qs = qi + q2 + Q, where Aw < ujj [j = 1, 2, 3). 
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Finally, taking into account that resonant mode interactions is a traditional 
subject of the nonlinear optics it is worth pointing out that the difference be- 
tween lattices and optical media comes from the anisotropy, which is an intrinsic 
property of the lattice. Also, lattices provide physical systems where almost all 
relations among parameters are available (including change of the signs of non- 
linear and dispersive terms). 
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